In-plane transport and enhanced thermoelectric performance in thin films of the 

topological insulators Bi2Te3 and Bi2Se3 



Pouyan Ghaemi,^'-^ Roger S. K. Mong,^ and J. E. Moore-^'^ 

^Department of Physics, University of California at Berkeley, Berkeley, CA 94720 
"^Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720 

Several small-bandgap semiconductors are now known to have protected metallic surface states 
as a consequence of the topology of the bulk electron wavefunctions. The known "topological insu- 
lators" with this behavior include the important thermoelectric materials Bi2Te3 and Bi2Se3, whose 
surfaces are observed in photoemission experiments to have an unusual electronic structure with a 
single Dirac cone [TJ[2].We study in-plane (i.e., horizontal) transport in thin films made of these 
materials. The surface states from top and bottom surfaces hybridize, and conventional diff'usive 
transport predicts that the tunable hybridization-induced band gap leads to increased thermoelec- 
tric performance at low temperatures. Beyond simple diffusive transport, the conductivity shows a 
crossover from the spin-orbit induced anti-localization at a single surface to ordinary localization. 

PACS numbers: 



Manipulation of phonon thermal conductivity and elec- 
tronic structure by quantum confinement [iSI have both 
been active areas of thermoelectric research, but the 
recent discovery that some of the best thermoelectrics 
are three-dimensional "topological insulators" (TIs) [H-E] 
suggests new directions for this field. For straightfor- 
ward reasons, similar materials are required for TI be- 
havior and for high thermoelectric figure of merit (FM) 
zT = ^-f-^, where S is thermopower (Seebeck) coeffi- 
cient, (7 is the electrical conductivity and k is the ther- 
mal conductivity. TIs require heavy elements in or- 
der to generate large spin-orbit coupling, which drives 
the formation of the topological surface state, and small 
bandgap so that the spin-orbit coupling can modify the 
band structure by "inverting" one band. 

Similarly, thermoelectric semiconductors typically con- 
sist of heavy elements in order to obtain low phonon ther- 
mal conductivity and have a small bandgap (of order 5-10 
times ksT, where T is the operating temperature) in or- 
der to obtain a large electronic power factor aS"^. An 
example is Bi2Te3, which has one of the highest known 
bulk thermopower FMs[7] and for decades has been used 
for thermoelectric devices [S]. 

Many approaches have already been tried to improve 
thermoelectric performance of Bi2Te3, e.g., tuning car- 
rier concentration (thereby increasing S), decreasing the 
thermal conductivity through alloying[9l [10], inducing 
crystal defects or nanostructures[lT], and making thin 
films to control transport of phonons and electrons [12]. 
Recent photoemission experiments reveal that Bi2Te3 is 
a TI with a single Dirac cone on the surface [2], consistent 
with electronic structure predictions [T^. It is natural to 
ask whether this new surface state could lead to improved 
thermoelectrics; this question is a major motivation for 
our study. 

We start by deriving the effective surface Hamiltonian 
on symmetry grounds, then use the resulting band struc- 
ture to compute the thermoelectric properties using stan- 



dard diffusive transport. Because of the large bulk gap 
(Agap ^ 165 meV), the chemical potential could be tuned 
to be in the bulk band gap]^, and in this regime we can 
study the transport properties of the surface states inde- 
pendently. In closing, anti-localization effects in a thin 
film resulting from the spin-orbit-induced Berry phase 
are discussed. 

The physical system to be studied here is a thin film 
of Bi2Te3. If the film is thin enough the surface states on 
both sides hybridize and open a gap |14lfl6] . Although 
the film thickness required to open an observable gap 
for TI surface states (^ 1-10 nm) is small, it is accessi- 
ble with current growth techniques and very recently the 
hybridization gap has been observed by in-situ photoe- 
mission on thin films of Bi2Se3. |17) 

We obtain the effective Hamiltonian of a thin film in 
the basis: |ut)i Wi), \dt), \d\.). The basis states are 
denoted by spin and the surface (top(u)/bottom((i)) at 
which they are localized. The general effective Hamilto- 
nian takes the block form: 

'^o' - ( vZ "hI ) 

Hu and are hermitian, while {VduY — Vud- The di- 
agonal block terms and Hd are the effective Hamil- 
tonians in the semi-infinite slab. The off-diagonal terms 
Vdu, Vud describe the coupling of the two surfaces. 

The parity (inversion) symmetry in the Bi2Te3 struc- 
ture puts constraints on the form of the Hamiltonian in 
(Eq[T]) that are strongest if the top and bottom surfaces 
are related by parity. Parity (H) exchanges the top and 
bottom surfaces, changes the sign of the momenta p, but 
has no effect on spin: H_ff(p)n^^ = iJ(— p). Time re- 
versal (7") exchanges the spins and flip the momenta: 
TH{p)T~'^ = H{-p). Explicitly H and T are 

n = (i M T^-^{'^\y)lc, 

where /C is the complex conjugation operator. 
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The combination of parity and time reversal implies 
nTi?(p)(nT)~^ = H{p), which restricts the form of the 
effective Hamiltonian via a^H^ay — a^H^ay — Hd and 
(jyVj^ay = Vdu- It useful note the identities ay{a')'^ay = 
— cr' for i = x,y,z and ayj^ay — I. If we decompose the 
matrix i7„ = cl + a. ■ cr, then (jyH^ay flips the sign of 
the Pauli matrices components a ct but leave the identity 
component cl: Hd = cl — a. ■ a. Since is hermitian, 
c is real and a is a real vector. A similar argument with 
Vdu = i!^fl + ■ <T will show that the Pauli matrices 
components must vanish, and that Vdu = Ay/ is a mul- 
tiple of identity, A/ is a complex number. 

''<p)^(itr-a('p'r.i)-<'" 

Without loss of generality, we can find a suitable gauge 
transformation such that A/ is real; this gauge transfor- 
mation is allowed by the phase ambiguity in the parity 
operator between the top and bottom states. 

Using parity to relate the parameters a, A/,c at op- 
posite momenta, there are further constraints: a(— p) = 
-a(p), A/(-p) = A/(p)* and c(-p) = c(p). Near the 
Dirac cone, we expand the Hamiltonian to the first order 
in p. As A^ and c are even in p, they can be treated as 
constant, while a is odd in p and we only keep the linear 
term, of the form VdP • c for Bi2Te3 (isotropic velocity). 

The resulting Hamiltonian is _ff = ^ yf' ^ , where 

Hu{Hd) corresponds to the upper (lower) surface of the 
film when the two surface wavefunctions have no overlap 
(thick film), and Vdu corresponds to the hybridization 
between the two surfaces. From first-principles electronic 
structure calculations [13j one can find parameters for H 
and the band structure of the surface modes: 

= ±^(VZ5P)2+A2, (4) 

Vd ~ 4.05 X 10^ to/s (in thin films is between 3.4 x 10^ 
and 4.6 x 10^ m/s [T7]) indicates the Dirac velocity and 
p is the momentum with respect to Dirac point. 

Now using the electron dispersion in Eq|4]we obtain the 
in-plane transport coefficients for the surface states [T5]: 

a ^ r de%^J:{e) (5) 

TaS = e p de^^^^^{e)[e^Ef) (6) 

T^o - f ds^^ne) {e - Eff (7) 

where as before a is conductivity and S Seebeck co- 
efficient, but kq is the electronic thermal conductivity 
at constant electrochemical potential, which is related to 
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FIG. 1: FM of surface states alone as a function of chemical 
potential and temperature. At low temperature the increase 
in the attainable FM becomes dramatic, but maximizing this 
increase requires good chemical potential control. 

the thermal conductivity for zero electric current through 
Kg = Ko — S'^aT. E (e) = N{e)T{e)iyxi£)^ is the conduc- 
tivity density and /(e) is the Fermi-Dirac distribution 
function. Experiments on bulk Bi2Te3 and its topologi- 
cally protected surface states allow estimation of the hy- 
bridization gap [2. .19j which we assume ^ 25meV. This 
could be achieved with recently created films of thickness 
a few nanometers J4l [151 Hi] • In Figure [I] we show the 
FM (zT) for the surface states alone. 

It can be seen that at low temperatures (i.e. 150 K and 
below) FM from surface states is large compared with 
current low temperature thermoelectric materials (e.g. 
Na^rCoOsHni and CsBi4Te6[lI]). In the thin film, trans- 
port properties of the surface states should be considered 
in parallel with the bulk of the thin film. Because of its 
thermoelectric performance {zT ~ 0.6) at room temper- 
ature, the bulk properties of Bi2Te3 have been widely 
studied[Hl[TS]. The FM for the whole film is 

{(Tf + Gb d) {kj + Kb d) 

f and b stand for surface and bulk properties and d de- 
notes the film thickness. Using the bulk properties [HI [53] 
for the temperature range where surface states have high 
FM, we can estimate the total thermoelectric FM for the 
thin film (but note that using a single value for the bulk 
contribution at a given temperature, since the transport 
properties even in bulk depend strongly on Fermi level). 
Figure 2 shows zT versus Fermi level for 1 nm thick films 
at 175 K, 150 K, 100 K and 50 K. The straight line in 
each plot indicates zT for CsBi2iTefi[2T]. 

Thermoelectricity from ballistic gapless ID edge states 
in a 2D quantum spin Hall (QSH) nanoribbon was dis- 
cussed recently by Takahashi and Murakami ^24j. Thin 
films of some 3D TIs are in QSH states [TS] , so the edge 
states could contribute also in a narrow film. However, 
in the QSH proposal the chemical potential is in the bulk 
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FIG. 2: FM for the thin film including bulk contributions 
at 50, fOO, f50, and 175 Kelvin. The straight line in each 
figure corresponds to the best known FM at the same tem- 
perature. Dashed line indicates the FM for the surface states 
alone and Dashed-Dotted hue indicates FM for bulk Bi2Te3. 
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bands; here it is in the gap. 

As is evident from Figure [2] at temperatures below 150 
K, which are important for several Peltier cooling appli- 
cations, the thermoelectric performance of the TI thin 
film is significantly enhanced because of the high FM of 
the protected edge states. At low temperature the bulk 
contribution is smaller than the surface contribution so 
that the unknown chemical potential dependence of bulk 
properties is not too significant. Crucially, the gap in the 
hybridized surface mode band structure can be controlled 
by tuning the thickness of the film to get high FM in a 
specific temperature range. The geometry of thin films 
is also very effective in reduction of phonon thermal con- 
ductivity, so there will be even larger enhancement for 
the TI thin films. The same approach applies as well to 
Bi2Se3 and Bi2Se2,Te3_a; alloys. 

The conductivity in EqjS] is calculated using Boltz- 
mann formalism [25] i.e. classical trajectories of elec- 
trons. Quantum mechanical corrections appear as the 
interference of quantum phases of different trajectories. 
In a simple picture, the amplitude for an electron mov- 
ing between two positions (a and b) is given by |Aa6| — 
I ^„ A„e*'^" I where n labels different classical paths be- 
tween a and b and (/)„ represents the corresponding quan- 
tum phase pS]. This phase is usually random and aver- 
ages to zero (i.e. |Aobp = J2n 1^"!^ )' which leads to 
the Boltzmann conductivity as in Eqj5] On the other 
hand, for closed trajectories where the electron returns 
to a region of size comparable to its quantum wavelength, 
quantum interference plays an important role. For dif- 
fusive motion, as in a metal, there is a finite probability 
for closed paths (fig. [3|. Such closed loops in ordinary 



metals lead to reduction of the conductivity (i.e. localiza- 
tion) [25]. The fact that such paths can be taken in two 
opposite directions means that the relative phases do not 
average to zero; the corresponding quantum correction is 
^ = — cos ((/)+ — 4>~) where (jf^ and 4>~ corresponds to 
phases for two opposite trajectories. 

Quantum corrections affect the conductivity of mate- 
rials strongly at mesoscopic length scales [27]. Here we 
show that the chiral nature of states on the Fermi sur- 
face (FS) can increase the mean-free path even when 
hybridization is present. At each momentum, the ef- 
fective surface Hamiltonian (Eqjsj) has two degenerate 
eigenstates: 
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where for simplicity = 1, p and B-p are the polar 
coordinates of p such that (pa;,py) =p (cos ^p, sin0p).On 
the FS -p is fixed so states are labeled by (?p's. 

Although at each momentum p states jop) and |6p) are 
orthonormal, they arc not at different momenta: 
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(6p|ap') is small for small angle scatterings and small A/. 
Any superposition of these two states (i.e. \ipp) = Q:|ap) + 

/3|6p)) is also an eigenstate with energy Ep = y^A^ + p'^. 

Smooth impurities on the surface scatter different 
states on the FS into each other (i.e. they change 9p) 
and lead to diffusive motion of quasiparticles and classi- 
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cal Boltzmann conductivity [25] a = ^^^^ which we used 
in our calculations. Quantum effects could be cast as 
quantum phases that quasiparticles gain while moving 
along different classical trajectories[2S]. Quantum cor- 
rections are due to interference of these different phases. 
However the phase difference between different trajecto- 
ries are generally random and on average give no con- 
tribution to the conductivity [25], with the exception of 
self-intersecting trajectories. Such trajectories could be 
taken in two different direction. The relative phase is not 
random and leads to quantum interferences. 

As mentioned above a general state on the FS is 
IV'p) = ci\ap) + P\bp) . But not all of these states can form 
a self intersecting trajectory. When Aj = (i.e. decou- 
pled surfaces) clearly the states \ap) is only affected by 
impurities in bottom surface and |6p) is affected by impu- 
rities in top surface. When impurities in top and bottom 




FIG. 4: Quantum corrections to the conductivity. Comparing 
with Fig. [2] we see a finite range over which antilocahzation 
coexists with enhanced zT predicted from ordinary diffusive 
transport (Eq. 5|6|7 1. 



FIG. 3: Back-scattering process: a) Two trajectories (A and 
B) leading to back scattering in real space traveled in opposite 
directions, b) A and B trajectories in momentum space. They 
add up to rotation of momentum by 2ti. 

surfaces are independently distributed, one can consider 
state I Op) or |6p) move over a self- intersecting trajec- 
tory independently, but not coherently together. In other 
words a mixed state of the form ji/ip) — a\a-p) + P\bp) 
with a ^ and ^ cannot go over a classical self- 
intersecting trajectory, since the direction of motion [Op) 
in I Op) and |5p) changes independently with scatterings in 
bottom and top surfaces respectively. As a result quan- 
tum interference effects are important for pure \a-p) and 
pure |6p) states not the mixed states. 

When Ay is non-zero but small, jop) and |6p) are 
mainly (but not entirely) localized in bottom and top sur- 
faces. An important feature of jop) and |6p) is that even 
when Ay ^ they only depend on direction of momen- 
tum (^p) in bottom and top surfaces respectively. Even 
though they have non-zero components on opposite sur- 
faces, they depend on the direction of momentum only 
in one surface. So as in the case with Ay = 0, when the 
impurities are randomly distributed on top and bottom 
surfaces, \a-p) and |&p) can move over self-intersecting 
trajectories independently. However, in a mixed state, 
\ipp) — ajflp) + (i\hp), the directions 9^ change indepen- 
dently in two surfaces and do not simultaneously form 
closed loops (see fig. [S]). As a result, in order to track 
the quantum interference effects we should consider jcp) 
and |6p) independently. The quantum phase of each tra- 
jectory is given by the associated Berry phase[28j: 




The quantum amplitude of the closed trajectory corre- 
sponds to the sum of the of amplitudes when the self- 
intersecting path is taken in opposite directions; i.e. 



I (e'"^! -I- 6'"^^) /y2|2 = 1 + cos - (f>2) where 0i and 
(/>2 correspond to the Berry phase of the passing over the 
closed loop in opposite direction. As can be seen in figure 
|3]((/)i — corresponds to the phase associated with ro- 
tation of 9p by 2n. The corresponding Berry phase could 

be calculated: 0i - (/.2 = rfi (*(pOlil*(pt)) = 

i^j^^Tr = ± ^"^2^^ TT for = I a) or 1 6). The two limits 

p ^ (A0 0) and p > Aj/Vd (A^ ±n) lead to 
localization and anti-localization behavior respectively. 
With increasing the chemical potential, localization be- 
havior changes into anti-localization, which is expected 
for a single TI surface. Comparing fig. |4]with the Seebeck 
results in figure [T] we see that in some range of chem- 
ical potential where the films have high thermoelectric 
performance, the anti-localization effect is also present. 
This leads to another advantage of the TI thin films: the 
quantum corrections increase the conductivity. 

We found that simply creating a nanometer-scale thin 
film of Bi2Te3 generates a hybridization- induced bandgap 
of the unconventional surface states that can be tuned by 
film thickness. Combined with some residual robustness 
against impurity scattering, this leads to an increased 
low-temperature FM of the film. The same technique 
can equally well be applied to Bi2Se3 or alloys in this 
class. While we have concentrated on the specific case of 
a thin film (or single superlattice layer) in this paper, the 
topological surface state of these materials is an impor- 
tant physical feature that will also affect thermoelectric 
transport in other nanoscale geometries. 
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